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ABSTRACT 

We have carried out 21 cm radio continuum, H76a radio recombination line, and 
various (^^CO, ^^CO, CS, & C^^S) molecular line observations of the W31 complex. 
Our radio continuum data show that W31 is composed of two extended HII regions, 
G10.2— 0.3 and G10.3— 0.1, each of which comprises an ultracompact HII region, two or 
more compact components, and diffuse envelope. The W31 cloud appears as an incom- 
plete shell on the whole and consists of southern spherical and northern flat components, 
which are associated with G10.2— 0.3 and G10.3— 0.1, respectively. For an assumed dis- 
tance of 6 kpc, the molecular cloud has a size of 48 pc and a mass of 6.2x10^ Mq- 
The IR luminosity-to-mass ratio and the star formation efficiency are derived to be 
9 Lq/Mq and 3%, respectively. These estimates are greater than average values of 
the inner Galactic plane. We detect two large (16 and 11 pc) and massive (2.1x10^ 
and 8.2x10^ -^0) CS-emitting regions in the northern and southern cloud components. 
The large amount (48% in mass and 16% in area) of dense gas may suggest that the 
W31 cloud has ability to form rich stellar clusters and that star formation has only 
recently begun. The extended envelopes of both G10.2— 0.3 and G10.3— 0.1 are likely 
to be results of the champagne flows, based on the distributions of ionized and molec- 
ular gas and the velocity gradient of II76a line emission. According to the champagne 
model, the dynamical ages of the two HII regions would be (4— 12) x 10^ yr. We find 
strong evidence of bipolar molecular outflows associated with the two ultracompact HII 
regions. In the vicinity of the ultracompact and compact HII regions in G10.3— 0.1, the 
^^CO J=2— 1/J=1— intensity ratio is high (1.4) and a small but prominent molecular 
gas hollow exists. Together these observations strongly indicate that the HII regions 
and their ionizing stars are interacting with the molecular cloud. Therefore, it is most 
likely that recently formed massive stars are actively disrupting their parental molecular 
cloud in the W31 complex. 
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1. Introduction 

Stars form in the dense cores of molecular clouds and eventually emerge as visible objects. 
Thus the molecular gas and dust surrounding newly formed stars should be disrupted. Massive 
stars are expected to play an important role in dispersing the ambient materials, since they emit 
strong radiation, stellar winds, and outflows even at their earliest evolutionary stage. Moreover, 
massive stars could dramatically impact on the surrounding environments through HI I regions and 
supernova explosions. However, direct observational evidence of this conjecture is rarely seen. 

W31 is one of the brightest HII region/molecular cloud complex in the inner Galaxy. Kalberla, 
Goss, & Wilson (1982) suggested on the basis of their sensitive HI absorption line observations that 
W31 is associated with the expanding 3-kpc arm with a peculiar velocity of ~30 km s~^ and so lies 
about 6 kpc from the Sun (see also Wilson 1974). We adopt this distance throughout the present 
paper, although there are some reports that propose different distances, e.g., 14.5ibl.4 kpc (Gorbel 
et al. 1997) and 3.4±0.3 kpc (Blum, Damineli, & Conti 2001). W31 appears as two prominent 
extended HII regions, G10.2— 0.3 and G10.3— 0.1, in the radio continuum maps with low (>1') 
angular resolutions (e.g., Reich, Reich, & Fiirst 1990). High-resolution (<10") radio continuum 
observations show that the central regions of both HII regions form elongated ionization ridges 
(Woodward, Heifer, & Pipher 1984; Ghosh et al. 1989). Two ultracompact (UC) HII regions, 
G10.15— 0.34 and G10.30— 0.15, were found to be located at the peaks of the ionization ridges, 
respectively (Wood & Churchwell 1989; Kim & Koo 2001). Several other massive young stellar 
objects were also identified in G10.2— 0.3 by near-infrared observations (Blum et al. 2001). These 
suggest that massive star formation is actively ongoing in the W31 complex. On the other hand, 
a Crab-like supernova remnant (SNR) GlO.0-0.3 is ~10' southwest of GlO.15-0.34. The SNR is 
likely to be associated with a soft 7-ray repeater, SGR 1806—20 (Kulkarni et al. 1994). Accordingly, 
the W31 complex seems to be one of the best laboratories for studying the formation of massive 
stars and the interaction of massive stars with their parental molecular cloud. 

In this paper, we present radio continuum, radio recombination line (RRL), and various molec- 
ular line observations of the W31 complex. Our observations aim to investigate the physical char- 
acteristics of molecular clouds that are actively forming massive stars, and to explore how massive 
stars and their natal molecular clouds interact with each other. The observations are described 
in § 2 and the results are presented in § 3. In § 4 we discuss the physical properties of the W31 
molecular cloud and the interaction between ionized and molecular gas in the region. Our main 
results are summarized in the last section. 



2. Observations 
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2.1. Radio Continuum and Recombination Lines 

The radio continuum and RRL data are extracted from the survey of 16 UC HII regions with 
extended envelopes. The observational details were presented by Kim & Koo (2001). In summary, 
radio continuum observations were made at 21 cm with the VLA of the NRAOQ in 1995 February. 
The array was in the DnC hybrid configuration and the observations were sensitive to structures 
up to 15'. The final map resulted from two separate pointings of the array. The phase centers were 
set at GlO.15-0.34 and GlO.30-0.15, i.e., (0,5)1950= (18^06°"22?5, -20°20'05") and (18^05""57?9, 
-20°06'26") or (a,<5)2000= (18^09'"21?0, -20°19'31") and (18^08'"56?1, -20°05'53"). The data 
were calibrated, imaged, and mosaicked using the AIPS. The resulting synthesized half-power 
beamwidth was 37" x 25" and the rms noise level of final image was 1.7 mJy beam~^. 

H76a (14.68999 GHz) RRL observations were conducted using the 140 foot telescope of the 
NRAO in 1997 February and June. The telescope has a FWHM of about 2' and a main beam effi- 
ciency of 0.58 at the observing frequency. Both circular polarizations were observed simultaneously 
using two 1024 channel autocorrelators with 40 MHz bandwidth each, yielding a velocity resolution 
of 1.59 km after Hanning smoothing. 

2.2. Molecular Lines 

We have carried out ^^CO, ^^CO, CS, and C^^S line observations of the W31 complex. A 
summary of the observations is given in Table 1. ^^CO J=l— and ^^CO J=2— 1 line observations 
were made with the 12 m telescope of the NRAO at Kitt Peak using the on-the-fly mapping 
technique. We used 256 channel filterbanks with 64 MHz and 128 MHz bandwidths as the backend, 
respectively. A squared area of 30' x 40' centered at (0,^)1950= (18''06"'09?1, -20°09'23") was 
mapped at 20" spacing in the ^^00 J=l-0 line, while a 5' X 5' region centered at G10.30— 0.15 was 
mapped at 10" sampling in the ^^CO J=2— 1 line. We used positions with some emission as reference 
positions, (/, b) = (10?15, — 1?50) and (10?05, -|-0?03), because we could not find any positions devoid 
of ^^CO and/or ^^CO emission near the source. We obtained a sensitive spectrum of each reference 
position and then added it to the position-switched spectra. The system temperatures were typically 
300 K at 110 GHz and 400 K at 230 GHz. We have converted the observed temperatures (T^) 
to the main-beam brightness temperature (Tb) using the corrected main-beam efficiencies (v^h) 
provided by the NRAO. 

CS J=2— 1 line observations were undertaken using the 14 m telescope of the Taeduk Radio 
Astronomy Observatory (TRAO). We mapped two regions of 11' x 14' and 12' x 10', respectively, 
centered at GlO.15-0.34 and GlO.30-0.15 in full-beam (60") spacing. The two UC HII regions 
were observed in the CS J=3— 2 and C^'^S J=2— 1 and J=3— 2 lines as well. We used SIS mixer 



^ The National Radio Astronomy Observatory is operated by Associated Universities, Inc., under cooperative 
agreement with the National Science Foundation 
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receivers with 256 channel filterbanks of 64 MHz bandwidth. The spectra were obtained by position 

switching to {I, b) = (10?15, 0?50), which was checked to be free from appreciable (T^ < 0.1 K) CS 
emission. The system temperature varied in the range 300—500 K at 98 GHz and 600—800 K at 
147 GHz. We also made ^^CO J=l— line observations of the region mapped in ^^CO 3=2—1 line 
emission using a 256 channel filterbank with 256 MHz bandwidth. The system temperature was 
about 1000 K. The line intensity was obtained on the scale from the chopper-wheel calibration. 

3. Results 

3.1. Radio Continuum and Recombination Line Results 

Figure 1 is our 21 cm radio continuum image of the W31 complex. There are two distinct 
extended HII regions, G10.2— 0.3 and G10.3— 0.1, in the field. G10.2— 0.3 has several compact 
components, almost all of which are embedded in low-level extended envelope of 10.'9 x 6.'7 (or 
19.0 pc X 11.1 pc). The central region is elongated in the east- west direction on the whole and is 
composed of two strong compact components. The UC HII region GIO. 15—0.34 lies at the peak 
of the western one. High-resolution (~5") VLA continuum observations at 5 GHz suggest that 
the central region comprises about 20 dense ionized clumps (Woodward ct al. 1984). The surface 
brightness declines sharply to the west, while it decreases slowly to the east. Figure 2 shows H76a 
spectra observed along 5(1950) ~ — 20°20'05". The center velocity decreases with increasing right 
ascension. The velocity difference is about 10 km s"^ between the central region and the eastern 
protuberance. 

In GIO. 3— 0.1, there are one isolated compact component and two strong ones surrounded 
by diffuse emission extending over 12.'8 x 4.'6 (or 22.3 pc x 8.0 pc). The two central compact 
components form an ionization ridge elongated in the northeast-southwest direction. The UC HII 
region GIO. 30— 0.15 is located at the peak of the western one. The diffuse envelope extends straight 
in the northwest-southeast direction, which is perpendicular to the central ridge. There is no 
significant velocity difference between the central region and the envelope. According to Kim & 
Koo (2001), the UC HII regions and their associated compact components are likely to be excited 
by the same ionizing sources in both GIO. 2— 0.3 and G10.3— 0.1, while the individual compact 
components may be ionized by separate sources. 

3.2. ^^CO J=l-0 Line Results 

Figure 3 exhibits ^^CO J=l-0 line profiles observed towards GlO.15-0.34 and GlO.30-0.15. 
There are three or more velocity components in the lines of sight. In each spectrum, the strongest 
component is associated with the UC HII region. The center velocities of ^^CO gas are shifted 
by 5—8 km s~^ from those of the ionized gas, marked by the vertical dotted lines. There seem 
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to be high-velocity wings in the ^^CO hnes of both UC HII regions, even though the features are 
significantly confused by multiple velocity components. Shepherd &: Churchwell (1996) have also 
found high-velocity wings in the ^^CO J=l— line of G10.30— 0.15, but could not identify them in 
GIO. 15—0.34 due to a complex blend of lines. We will discuss in more detail the high- velocity gas 
using our CS spectra in the next section. 

Figure 4 is ^^CO line integrated intensity map, obtained by integrating over the velocity range 
between ulsr = and 22 km s~^. Large crosses indicate the positions of UC HII regions, while small 
crosses represent the peak positions of the strong compact components without UC HII regions (see 
Fig. 1). The molecular cloud is very clumpy. There is a large hollow over the central region, which 
makes the cloud appear as an incomplete shell being open to the northeast. A possible explanation 
for this morphology is that the cloud has been produced by energetic events near the center, such as 
strong stellar winds and/or supernova explosions of massive stars, although we have not found any 
massive star or nonthermal radio continuum emission in the hollow region. As will be discussed in 
§ 4.2, the dynamical age of G10.2— 0.3 is similar to that of G10.3— 0.1, suggesting that star formation 
has begun almost simultaneously at the southern spherical and northern flat components. This may 
be in favor of the above interpretation. It seems that the northern component has been substantially 
dispersed by massive stars, since a small but conspicuous hollow exists near the UC and compact 
HII regions, (q;,5)i95o ~ (18'^05™58**, — 20°05'45") (see § 4.2). In the southern component, on the 
contrary, the dense ^^CO cores associated with the UC and compact HII regions are still strong. 
It is worthwhile to note that the integrated intensity drops steeply on the eastern boundary of the 
southern component, in contrast to the radio continuum intensity distribution. 

Figure 5 displays the channel maps of the molecular cloud. ^^CO emission appears mainly in 
the velocity range between 2 and 20 km s"-*^. The channel maps show that the northern and southern 
components are certainly associated, and reveal more clearly incomplete shell-like morphology of the 
cloud. The molecular gas shell remains nearly constant in size as the velocity varies, suggesting that 
it does not dynamically expand or contract. There are small velocity variations in both northern and 
southern components. As the velocity decreases, strong ^^CO emission moves from the western part 
to the eastern part in the northern component, while it moves from the northeastern part through 
the northern part to the southern part in the southern component. The ^^CO core associated with 
GIO. 15— 0.34 (--^8 km s~^) is slightly blueshiftcd from the main shell-like structure (~12 km s~^). 
The northern flat component seems to be nearly broken through at the position of the distinctive 
hollow, which persists over a velocity interval of 10—15 km s^^. 

Assuming that the molecular gas is in local thermodynamical equilibrium (LTE), ^^CO column 
density, iV(^^CO), can be computed from the formula (see, e.g., Dickman 1978) 

m) = 2.42xl0"^-^ll^|^ (cm-), (1) 
where Tgx is the excitation temperature in K, r is the optical depth, v is the velocity in km s~^. The 
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mass of molecular cloud, Mlte, is measured to be 6.20x10^ M©, provided that -"^^CO abundance 
relative to H2 is 2xlO~^ (Dickman 1978). In this calculation we assumed that the cloud boundary 
is / T^dv ~ 10 K km s~^ and that T^x is 30 K in the central regions of the northern and southern 
components and 15 K in the outer regions, based on our ^^CO line data and the Massachusetts- 
Stony Brook ^^CO Galactic Plane Survey (Sanders et al. 1986). Our mass estimate adopted a 
mean molecular weight of 2.3 to account for the contribution of helium. Table 2 summarizes the 
physical parameters of the W31 molecular cloud. 

3.3. CS Line results 

Figure 6 exhibits our CS and C^^S line profiles obtained towards GlO.15-0.34 and GlO.30-0.15. 
The line parameters of these spectra are presented in Table 3. There is a single velocity component 
along the line of sight and it is associated with each UC HII region. The center velocity of CS gas is 
in good agreement with that of ^^CO gas. The CS 3=2—1 and J=3— 2 lines of both UC HII regions 
show distinct high-velocity wings. This implies that dense gas exists in high-velocity molecular 
components. There are a few tens of massive star-forming regions where high-velocity wings have 
been observed in the CS lines. For example, Plume, Jaffe, & Evans (1992) detected CS J=7— 6 line 
emission in 104 massive star-forming regions, selected originally by the presence of an II2O maser, 
and found high- velocity wings in more than 18 out of them. The full widths of the CS J=2— 1 
lines are about 50 km s^^ at typical rms noise level for GIO. 15—0.34 and about 30 km s^^ for 
G10.30— 0.15. These values are similar to those expected from the ^^CO lines, respectively. The 
CS J=3— 2 line of GIG. 30— 0.15 has very similar shape and full width as the CS J=2— 1 line, while 
the CS J=3— 2 line of GIO. 15— 0.34 has a different asymmetric shape and a smaller (20 km s""*") full 
width than the CS J=2-l line. The fuh width of the CS J=2-l line of GIO. 15-0.34 is comparable 
to that of G5. 89— 0.39, which drives one of the most energetic outflows in the Galaxy (Harvey &: 
Forveille 1988). On the other hand, we can not find high- velocity wings in the CS lines observed 
in the surrounding regions of both UC HII regions. It seems mainly because the extents of CS 
outflows are ;^1', taking into account that CO outflows observed in massive star-forming regions 
usually have linear sizes of ^2 pc (Ridge & Moore 2001; Beuther et al. 2002). So high-resolution 
CS line observations are required to confirm whether the high-velocity gas is indeed due to bipolar 
molecular outfiows, and to explore the physical and dynamical properties of the outflows. 

Assuming that the excitation temperature and the beam fllling factor are identical in the CS 
and C'^^S lines in each transition, wc determine the optical depth of C^^S line emission, Tp, from 
the peak brightness temperature ratio using the following formula 

Tb(C34s)/rb(CS) = [1 - exp(-rp)]/[l - exp{-TpR)]. (2) 

The estimated values are listed in Table 3. Here the CS to C^^S abundance ratio, R, was assumed 
to be equal to the terrestrial value (22.5). 
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Figure 7 shows integrated CS 3=2—1 line intensity maps. The CS gas distribution agrees well 
with the ^^CO gas distribution in the dense regions of the northern and southern components. We 
can see the prominent hollow near G10.30— 0.15 in the CS gas distribution as well. The UC HII 
regions correspond to the peaks of CS cores, while the compact components without UC HII 
regions do not. This seems to be consistent with the result of Kim & Koo (2001) that the compact 
components with UC HII regions are smaller and denser than those without UC HII regions. 
Based on the observation they suggested that the former are in an earlier evolutionary stage than 
the latter. We have also found the same trend in several other UC HII regions with extended 
envelopes that have two or more compact components. The results will be presented in a separate 
paper (Kim & Koo 2002). 

If CS line emission originates from regions in LTE and is optically thin, the CS column density, 
A'"(CS), is given as 



where Tr is the rotation temperature in K, Tb is the brightness temperature in K, and v is the 
velocity in km s"^. If we take a fractional CS abundance to H2 of 1x10"^ (Linke & Goldsmith 

1980; Frcrking et al. 1980), the total masses of CS-emitting regions are 2.13x10^ Mq for GlO.2-0.3 
and 0.82x10^ Mq for GlO.3-0.1. In these calculations Tr was assumed to be 30 K on the basis of 
the results of studies on dense cores in massive star-forming regions (e.g., Snell et al. 1984; Linke 
& Goldsmith 1980). Table 2 lists the physical parameters of the CS-emitting regions. 



We obtained ^^CO J=2— 1/J=1— integrated intensity ratio map for the squared area of ~5'x5' 
centered at GlO.30-0.15, after convolving ^^CO J=2-l data to the resolution (53") of ^^CO J=l-0 
data. The integrated velocity range was vlsr = 0—22 km s~^. Figure 8 shows the distribution of 
^^CO J=2— 1/J=1— ratio (grey scale) overlaid with the distribution of ^^CO J=2— 1 line integral 
(contours). The ratio peaks over the dense cores associated with the UC and compact HII regions 
and is low (^1.2) over the hollow region. The estimated values are between 1.0 and 2.0. The 
average value is 1.4, which is much greater than typical value (0.5—0.8) for molecular clouds in the 
Galactic disk (e.g., Sakamoto et al. 1997). Such a high (>1.0) ratio is observed towards the areas 
around HII regions in molecular clouds undergoing massive star formation, e.g., Orion A molecular 
cloud (Castets et al. 1990; Sakamoto et al. 1994), Sgr A molecular cloud (Oka et al. 1996), and 
W51B molecular cloud (Koo 1999). 

If a molecular cloud has a uniform temperature, these high ratios can be obtained when the 
molecular gas is warm, dense, and less opaque. For example, Sakamoto et al. (1994) showed on 
the basis of an one-zone large velocity gradient (LVG) analysis that ratios greater than unity are 




(3) 



3.4. 



CO J=l-0 and J=2-l Line Results 
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produced when the CO-emitting regions have rex;^40 K, n(H2);^lxlO^ cm~^, and T2i<5. In this 

model the ratio variation is mainly attributed to the density variation of unresolved clumps that 
constitute the cloud. However such a single-temperature model does not seem to be applicable 
to our case, since the ratio of ^^CO J=l— and ^"^CO J=l— line intensities, ^^/^^i?i_o, ranges 
between 1.3 and 6.7 with an average of 2.7, indicating that ^^CO line emission is optically thick 
(r2i»5). On the other hand, the high ratios can also be obtained when a dense molecular cloud 
is externally heated by strong far-ultraviolet (FUV) radiation. Gierens, Stutzki, & Winnewisser 
(1992) performed a detailed ^^CO and -^^CO line radiative transfer analysis for spherically sym- 
metric clumps irradiated by intense FUV radiation, and showed that the ratios > 1 observed in 
the Orion A molecular cloud are easily explained by the temperature and abundance variations 
within the photodissociation region layer of the clumps. If we adopt their results, we can de- 
rive <n(H2)> and <A''(H2)>/Ai; of the CO-emitting regions from the ratio of ^^CO J=2— 1 and 
J=l— line intensities, ^^i?2-i/i-05 ^^^^^Ri-o using their Figures 3 and 5. The ^^-^2-1/1-0'^ 
are 1.0—2.4 with an average of 1.6. The average values for ^"^^^^Ri^q and ^^^2-1/1-0 suggest 
that <iV(H2)>/Aw;^lxl0^2 cm'^ (km s"^)"^ and <n(H2)>~10^ cm"^. The ^^i22-i/i-o is lower 
(^1.4) over the hollow region than the surrounding regions, while the ^^^^^Ri-q is higher (^4). 
This implies, as expected, that <n(H2)> and <A''(H2)>/Ai; are both lower in the hollow region, 
i.e., <7V(H2)>/At;=(0.3-1.0)xl022 cm'^ (km s-i)-^ and <n(H2)>=(1.5-6.0)xl03 cm-^. 



4. Discussion 

4.1. Physical Characteristics of the W31 Molecular Cloud 

The W31 molecular cloud is classified as a giant molecular cloud (GMC) based on the physical 
parameters determined from the ^'^CO line observations. The CS-emitting regions of W31 are 
much more massive and larger than those (10^—10^ Mq in mass and 0.5—5 pc in size) of molecular 
clouds associated with Sharpless HII regions, which are excited by single late O or early B stars 
(Zinchenko et al. 1994; Carpenter, Snell, Sz Schlocrb 1995). Hence the fraction of dense (>10^ 
cm~^) molecular gas may be much higher in the W31 cloud than in the molecular clouds with 
Sharpless HII regions. For W31 as a whole, the CS-emitting regions contain 48% of the total mass 
and occupy 16% of the total area (see Table 2). In comparison, these values are, respectively, one 
order of magnitude greater than estimates (2% and <1%) for the Gemini OBI cloud complex with 
about 10 Sharpless HII regions, which has a similar total mass to the W31 cloud (Carpenter et al. 
1995). This seems to be consistent with the observation that more luminous IRAS sources tend 
to be associated with more massive CS cores in the Gemini complex, suggesting that more massive 
cores in general form more massive stars (Carpenter et al. 1995). Similar results have been found 
in the Orion B (L1630) molecular cloud (Lada 1992) and the Rosette molecular cloud (Phelps & 
Lada 1997). In these molecular clouds the embedded clusters are associated with the most massive 
dense cores. These observations indicate that both high gas density and high gas mass are required 
for the formation of stellar clusters. On the other hand, the CS-emitting regions of W31 are further 
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more massive and larger than those of some other GMCs undergoing active massive star formation, 
such as M17 (Snell et al. 1984) and Orion B (Lada et al. 1991). This difference is Ukely due 
to difference in evolution (destruction of dense cores by massive stars) as well as difference in the 
initial physical properties, since some newly formed massive stars have already emerged as visible 
objects in the regions (e.g., Stephenson Sz Hobbs 1961; Warren Sz Hesser 1978). 

The ratio of infrared (IR) luminosity to mass is a good measure of the massive star formation 
activity of molecular clouds. The ratio varies widely from cloud to cloud, and shows little correlation 
with cloud mass from 10^ to 10^ Mq or location in the Galaxy (Mooney &; Solomon 1988; Scoville 
&; Good 1989; Carpenter, Snell, &; Schloerb 1990). We estimate the ratio for the W31 molecular 
cloud using IRAS High RESolution processing (HIRES) images at 60 /xm and 100 /Ltm. After a 
smooth background emission is removed from each image, the total flux densities at 60 ^m and 
100 /xm, Feo and Fioo, are measured to be (7.4 ± 0.8) x 10^ Jy and (1.09 ± 0.10) x 10^ Jy within 
the could boundary, respectively. Here we took average intensities on the cloud boundary as the 
background emission levels. The total IR luminosity in the wavelength range 1—500 /xm, -Zjir, can 
be calculated from Fqq and Fioo in Jy using the following formula (Lonsdale et al. 1985; Lee et al. 
1996) 

LiR = 0.394 R{%,(3) [ Fioo + 2.58 Fqq ] (f {Lq), (4) 

where is the 60/100 /^m color temperature in K, /3 is the index in the emissivity law, Qabs ~ 
and d is the source distance in kpc. The R{T^,f3) is color correction factor that accounts for the 
flux radiated outside the 60 jiva and 100 /xm IRAS bands. If (3 is assumed to be 1 (Hildebrand 

1983), LiR is derived to be 5.8x10^ Lq, which gives Lir/Mlte=9.4 Lq/Mq. This estimate is much 
greater than the average ratio (2.8 Lq/Mq) of all molecular clouds in the inner Galactic plane, but 
is comparable to the median value (7 Lq/Mq) for GMCs with HIT regions (Scoville & Good 1989). 

We also derive star formation efficiency, SFE= M*/(M* + Afdoud); of the W31 cloud. The 
stellar mass was computed from the Lynman continuum photon flux measured by radio 
continuum observations of GIO. 2-0.3 and GlO.3-0.1, 2.7x10^° photons s'^ (Kim & Koo 2001). 
Here we adopt < M* > / < >= 5.7 x 10"^^ Mq (photons s~^)~^, which was determined by 
McKee & Williams (1997) using the initial mass function of Scalo (1986) and the parameters of OB 
stars of Vacca, Gar many, Sz Shull (1996). The SEE estimated for the entire cloud is 3%, which is 
somewhat greater than the Galactic median value of ~2% (Myers et al. 1986). On the other hand, 
the SFE's for the two CS-emitting regions are derived to be 6% and 4%, respectively (Table 3). 

4.2. Interaction between Ionized and Molecular Gas 

Figure 9 shows the ionized gas distribution (contours) superposed on the molecular gas distri- 
bution (color). For both GIO. 2— 0.3 and GIO. 3— 0.1, the morphology of extended envelope seems 
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to be largely determined by the distribution of the ambient molecular gas. In G10.2— 0.3, the radio 

continuum surface brightness decreases slowly to the cast while it drops rapidly to the west. This 
surface brightness behavior can be naturally understood if the central part of G10.2— 0.3 is situated 
within the molecular cloud and the eastern protuberance is a result of the champagne flow (Israel 
1978; Tenorio-Tagle 1979). A large (>10 km s~^) velocity gradient exists between the central and 
eastern parts (see Fig. 2), which is compatible with what would be expected from the champagne 
flow model. On the other hand, the extended envelope of G10.3— 0.1 has a bipolar morphology. 
Such a morphology can be created by the champagne flow in case where the HII region forms in 
a flat molecular cloud (Bodenheimer, Tenorio-Tagle, & Yorke 1979; Tenorio-Tagle, Yorke, Boden- 
heimer 1979). Indeed the distribution of molecular gas in the region is severely elongated in the 
direction orthogonal to the bipolar axis of G10.3— 0.1 (see also Fig. 7b). Therefore, the extended 
envelope of G10.3— 0.1 is also likely to be a consequence of the champagne flow. However, we have 
not observed any significant velocity difference between the central region and the envelope. This 
might be due to a small inclination. The velocity of a champagne flow is determined by the pressure 
contrast at the cloud boundary and can increase up to about 30 km s""*^ (Tenorio-Tagle 1979). Since 
the farthest boundary is ~12 pc (7') from the central region for both G10.2— 0.3 and G10.3— 0.1, 
their dynamical ages may be in the range (4—12) x 10^ yr. For G10.2— 0.3, our estimate matches 
with the result of Blum et al. (2001), who derived a mean age of ~lxlO^ yr for the embedded 
cluster associated with the eastern compact component in the central region, based on their near-IR 
photometry and spectroscopy. 

A striking feature is the distinctive molccTilar gas hollow in the northern flat component of 
the W31 cloud. The hollow is ~4 pc (2.'4) in size. It is located close to the UC and compact HII 
regions. The thickness of the flat molecular cloud is much smaller at the hole position than other 
parts. These observations strongly suggest that the hollow might have been produced due to the 
dispersion and destruction of molecular gas by the HII regions and their ionizing stars. This view is 
supported by the fact that the ^^CO J=2— 1/J=1— integrated intensity ratio exceed unity around 
the UC and compact HII regions, implying that the HII regions and their central stars are closely 
interacting with the surrounding molecular materials. It is well known that the surrounding gas 
and dust of newly formed massive stars should be eventually dispersed and destroyed by strong 
radiation, stellar winds, and HII regions. 

If the molecular gas column density of the hollow region had been equal to the average value 
of the surrounding regions, which have / T^dv ^ 25 K km s~^, the amount of molecular gas dis- 
persed by the massive stars would be (4ibl)xl0^ Mq. The quoted error represents variation in 
the measured mass depending on the assumed threshold value for the surrounding regions. Our 
estimate is not very different from the ionized gas mass of GIO. 3-0.1, 2x10^ Mq, determined from 
radio continuum observations (Kim & Koo 2001). The derived mass and dynamical age imply 
that the disruption rate of the molecular cloud is (3— 10) x 10^^ Mq yr^"*^, which is very similar 
to the results of two-dimensional numerical calculations for the champagne flows in various condi- 
tions, (5-10) X 10-^ Mq yi-^ (Yorke, Bodenheimer, & Tenorio-Tagle 1982). Whitworth (1979) also 
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performed a semi-analytic study on the efficiency with which massive stars destroy their parental 
molecular clouds and showed that the disruption rate is determined by the Lyman continuum pho- 
ton flux, the ambient density, and the elapse time. The Lyman continuum photon flux of G10.3— 0.1 
is 5.5x10'^^ photons s~^ (Kim Sz Koo 2001). The average density of the surrounding regions with 
/ T^dv ;^ 25 K km s"^ is ~lxlO^ cm ^. In that case the disruption rate is (5±l)xl0 ^ M© yr ^. 
This value is in good agreement with the rate estimated above, too. In addition to strong radiation, 
the bipolar outflow associated with G10.30— 0.15 might have played a role in the dispersion pro- 
cess. Therefore, it is most likely that the northern flat molecular cloud is being disrupted by newly 
formed massive stars. This may also be true for the southern spherical molecular cloud, although 
compelling observational evidence for molecular gas disruption have not been found except the 
champagne flow. 

5. Conclusions 

We have carried out 21 cm radio continuum, H76a RRL, ^^CO, ^=^C0, CS, and C^^'S line 
observations of the W31 HII region/molecular cloud complex. The physical parameters of the 
W31 cloud arc largely comparable to those of other GMCs that are undergoing active massive 
star formation. However, the CS-emitting regions are more massive and larger than those of some 
other GMCs with indicators of massive star formation. The large amount of dense (>10'^ cm^'^) 
gas implies that the W31 cloud has ability to form rich stellar clusters and suggests, together 
with the young dynamical ages (;^lxlO^ yr) of G10.2— 0.3 and G10.3— 0.1, that star formation has 
recently begun in the cloud. On the other hand, we have found high-velocity molecular gas towards 
GIO. 15—0.34 and G10.30— 0.15, which are likely due to bipolar outflows. The extended envelopes of 
both GIO. 2— 0.3 and GIO. 3— 0.1 may be results of the champagne flows, based on the distributions 
of ionized and molecular gas and the velocity gradient of RRL emission. In the vicinity of the UC 
and compact HII regions in GIO. 3— 0.1, a distinctive molecular gas hollow is present and the ^^CO 
J=2— 1/J=1— line intensity ratio exceeds unity. Together these observations support that the HII 
regions and their ionizing stars are strongly interacting with the surrounding molecular materials. 
Therefore, W31 is a good example showing not only formation of massive stars but also destruction 
of their parental molecular cloud by them. 

We are grateful to Bob Blum, John Carpenter, Youngung Lee, and Young Choi Minh for 
helpful comments, and Yong-Sun Park for giving us part of his observing time at the TRAO. We 
also thank the staff at the IPAC in Caltech for processing the IRAS HIRES maps for W31 region. 
This work has been supported by a bilateral KOSEF-Korea and CONACyT-Mexico agreement 
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Fig. 1. — 21 cm radio continuum image of W31 made with the VLA (DnC array). This image 
is corrected for the primary-beam attenuation. The synthesized beam is 37" x 25" at P.A.=76°. 
Contours are 10, 30, 70, 150, 300, 600, 1000, 1500, and 2300 mJy beam'^ Large crosses are the 
positions of GIO. 15—0.34 and G10.30— 0.15, while small crosses are the positions where H76q! line 
emission was observed. The SNR GIO. 0-0. 3 is marked by an eight-pointed star and the Galactic 
midplane is represented by the dotted line. A 5 pc linear scale bar is shown in the bottom left 
corner. 



_l J I L_ 



_l I I L_ 



-50 



50 



100 



VrsR (km s-0 



Fig. 2.— H76a line spectra observed along ^(1950) ~ -20°20'05" in GlO.2-0.3. Right ascension 
offsets with respect to GIG. 15—0.34 are given in second on the right side. The vertical dotted line 
marks the center velocity of H76a line emission observed toward GIG. 15— 0.34, i;lsr=16.4 km s~^. 
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Fig. 3.— i^CO J=l-0 line profiles obtained towards (a) GlO.15-0.34 and (b) GlO.30-0.15. In 
each panel, the vertical dotted line indicates the center velocity of 1176a line. 
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Fig. 4. — '^CO J=l— integrated intensity map. The integrated velocity range is from wlsr=0 to 
22 km s-^ Contour levels are 13, 20, 25, 30, 40, 50, 60, 70 and 80 K km s"^ and grey scale flux 
range is 10—86 K km s"^. The boxes show the areas covered by our CS J=2— 1 line observations. 
Large and small crosses represent the UC HIT regions and strong compact components without 
UC HIT regions, respectively. The SNR GIG. 0-0. 3 is marked by an eight-pointed star and the 
Galactic midplane is indicated by the dotted line. 
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Fig. 5. — Channel maps of CO J=l— line emission. The center velocity is given in km s at 
the upper-right corner in each panel. Contours are 1, 2, 3, 4, 6, and 8 K and grey scale flux range 
is 0.5—10.0 K. Symbols are the same as in Figure 4. 
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Fig. 6. — CS J=2— 1 and J=3— 2 and C^^S 3=2—1 and J=3— 2 line profiles observed towards (a) 
GlO.15-0.34 and (b) GlO.30-0.15. The center velocity of H76q; line is marked by the vertical 
dotted line in each panel. 
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Fig. 7.— CS J=2-l integrated intensity maps of (a) GlO.2-0.3 and (b) GlO.3-0.1. Contour 
levels are 1.5, 3, 4, 5, 6, 7.5, 9, 11, and 14 K km s"^ in (a), and 1.2, 2, 3, 4, 7, 12, and 18 K km s"^ 
in (b). Large and small crosses represent the UC HII regions and strong compact components 
without UC HII regions, respectively. 
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Fig. 8. — ^^CO J=2— 1/J=1— integrated intensity ratio map (grey scale) and ^^CO J=2— 1 
integrated intensity map (contours) of a '~5'x5' area around G10.30— 0.15. Grey scale flux range 
is 1.2—2.0 and contours are 80, 100, 120, 140, 170, and 240 K km s~^. Large cross represents 
G10.30— 0.15, while small cross indicates the strong compact components without UC HII region. 
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Fig. 9. — ^^CO integrated intensity map (color) and 21 cm radio continuum image (contours). 
The integrated velocity range is 0-22 km s'^. Contour levels are 10, 30, 70, 200, 500, 1000, and 
2200 mJy beam~^. The color wedge on the right shows the integrated intensity range. Crosses 
represent the positions of G10.15— 0.30 (upper) and G10.34— 0.15 (lower). 



Table 1. MOLECULAR LINES OBSERVED 



Frequency Resolution FWHM OFF Observation 



Telescope and Transition (GHz) A/ (kHz) Av (km s ) (") i^mb H, b) Date 



NRAO 12-m 

^=^CO J=l-0 110.20135 250 0.68 60 0.84'' (10?15, -1?50)'' 1997 June 

^^CO J=2-l 230.58300 500 0.62 27 0.51*^ (10.05, +0.30)'' 2000 March 

TRAO 14-m 

OS J=2-l 97.98095 250 0.77 60 0.48 (10.15, +0.50) 1997 March 

CS J=3-2 146.96904 250 0.51 46 0.39 (10.15, +0.50) 2000 March 

C^^S J=2-l 96.41298 250 0.77 60 0.48 (10.15, +0.50) 2000 March 

C^*S J=3-2 144.61711 250 0.51 46 0.39 (10.15, +0.50) 2000 March 

^^CO J=l-0 115.27120 1000 2.60 53 0.46 (10.05, +0.30)'' 2000 March 



'^Corrected main beam efficiency r;*^^, Th—T^/rj^i^. 
''Reference positions with some '^'^CO and/or ^"^CO emission. 



-25- 



Table 2. ^^CO J=l-0 AND CS J=2-l LINE OBSERVATIONAL RESULTS 







R 


Mlte 


n(H2) 




log iV^'^ 




SFE 


Transition 


Source 


(pc) 


(10* Me) 


(10^ cm"^) 


(km s-i) 


(s-^) 


(10* Mq) 


(%) 


i^CO J=l-0 


W31 


24.1 


62.0 


2.0 


8.2 


50.43 


1.54 


3 


CS J=2-l 


GlO.2-0.3 


8.0 


21.3 


19.2 


10.0 


50.33 


1.22 


6 




GlO.3-0.1 


5.6 


8.2 


21.7 


5.4 


49.74 


0.31 


4 



TWHM of the average spectrum of the entire mapping area. 
''Taken from Kim & Koo 2001. 
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Table 3. CS AND C^^^'S LINE PARAMETERS 



V-LSR Tl Av Tp 

Source and Transition (km s~^) (K) (km s~^) 



GlO.15-0.34 

CS J=2-l 9.6 1.31 8.8 

CS J=3-2 10.8 1.05 7.1 

C^^S J=2-l 10.4 0.32 3.1 0.28 

C'^^S J=3-2 10.3 0.50 2.6 0.65 

GlO.30-0.15 

CS J=2-l 13.5 1.90 8.8 

CS J=3-2 13.8 2.00 7.7 

C^^S J=2-l 15.1 0.40 7.0 0.24 

C^*S J=3-2 13.4 0.62 2.1 0.37 
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